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Abstract 

We investigate various aspects of a proposal by Aharony and Kutasov [1] for the gravity dual 
of an open Wilson line in the Sakai-Sugimoto model or its non-compact version. In particular, 
we use their proposal to determine the effect of finite temperature, as well as background 
electric and magnetic fields, on the chiral symmetry breaking order parameter. We also 
generalize their prescription to more complicated worldsheets and identify the operators 
dual to such worldsheets. 
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1 Introduction 



In [1] Aharony and Kutasov gave a prescription for computing correlators of open Wilson line 
(OWL) operators in the Sakai-Sugimoto model and its non-compact version. These OWLs 
are interesting since their vevs are order parameters for chiral symmetry breaking (%SB). 
Recall that the Sakai-Sugimoto model [2] is holographically dual to a gauge theory, known as 
holographic QCD, which shares many dynamical features with QCD- namely, confinement 
and xSB. The non-compact Sakai-Sugimoto model is holographically dual to a theory which 
we will call holographic Nambu-Jona-Lasinio (NJL) which, like the usual NJL model, has 
X;SB [3]. We review some details of these models below. 

In this paper, we use the prescription of [1] to probe how the x^B order parameters 
in holographic QCD and NJL depend on temperature, background electric and magnetic 
fields, and various parameters defining the OWL operator. Since confinement is not our 
main interest, we mostly focus on holographic NJL, where confinement is turned off. Most 
of our results are qualitatively unchanged if we consider, instead, holographic QCD. This is 
analogous to the case in field theory where, when it comes to the analysis of chiral symmetry 
breaking, the NJL model shows qualitatively similar behavior to QCD. 

We will first give a brief review of the string models for holographic QCD and NJL, and 
for the proposals for %SB order parameters in these models. The Sakai-Sugimoto model [2] 
is an intersecting brane model made of N c "color" D4-branes extended in the (x°x 1 x 2 x 3 x' 1 )- 
directions, with x 4 being a Scherk-Schwarz circle of radius R intersecting jV f D8-branes and 
Nf D8-branes at two (3 + l)-dimensional intersections. The "flavor" D8- and D8-branes are 
located at the antipodal points on the x 4 -circle. Imposing anti-periodic boundary conditions 
for fermions around the x 4 -circle leaves the U (N c ) gauge boson modes of the 4 — 4 strings 
massless, but gives mass to their fermionic and scalar modes. There are also massless Weyl 
fermion modes localized at the D4-D8 and D4-D8 intersections, denoted ip L and ip R , which 
come from the Ramond-Ramond sector of the 4-8 and 4-8 strings, and which transform as 
(N c ,N f ,l) and (N c ,l,N f ) of U(N C ) x U(N { ) x U(N { ), respectively. The U(N { ) x U(N { ) 
gauge symmetry of the flavor branes is the chiral symmetry for these Weyl fermions. 

The low energy theory on the color branes is a (4 + l)-dimensional U (N c ) gauge theory 
with a dimensionful 't Hooft coupling A5. The dynamics of the theory at the intersections 
are governed by the dimensionless effective 't Hooft coupling A 4 = A 5 /i? - At weak cou- 
pling (A 4 <C 1), the low energy effective theory contains QCD but is hard to analyze using 
gauge-gravity duality [4, 5, 6, 7], while at strong coupling (A4 ^> 1) it can be effectively 
analyzed this way; however, the theory is no longer QCD but is instead a related theory 
known as holographic QCD, with similar qualitative features. At strong coupling, and in 
the probe approximation JVf <C N c , one can safely consider the flavor branes in the near- 
horizon geometry of "color" D4-branes. The analysis of the flavor DBI action shows that 
the energetically-favored configuration of the flavor branes is one where they smoothly join 
at some radial point in the background geometry to form a U-shaped configuration where at 
the tip of the configuration the asymptotic U(N{) x U(N{) chiral symmetry is broken down 
to a single U(N{). By analyzing the near-horizon geometry of the color branes, one can also 
show that the model exhibits confinement and a mass gap for glueballs [8, 9]. 
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One can separate the scale of chiral symmetry breaking from that of confinement by 
allowing the flavor branes to be separated by an asymptotic distance of £q < ttRq, instead 
of being placed at the antipodal points of the x 4 -circle as in the Sakai-Sugimoto model [3]. 
One can completely turn off confinement by taking the Rq — > oo limit where the field theory 
dual is holographic NJL (which reduces to a nonlocal version of the NJL model at weak 
coupling). At strong coupling, chiral symmetry is broken via a smooth fusion of the flavor 
branes, providing a holographic model of %SB without confinement. 

These models can also be considered at finite temperature and chemical potential, where 
they exhibit similar behavior to QCD at finite temperature and chemical potential. The 
finite-temperature analysis of the Sakai-Sugimoto model [10] and its non-compact version 
[11] amounts to putting the probe flavor branes in the near-horizon geometry of iV c non- 
extremal D4-branes. The DBI action of the flavor branes indicates that at low temperatures 
(compared to i^ 1 ) the energetically-favored solution is that of smoothly-connected D8- and 
D8-branes which, like its zero-temperature counterpart, is a realization of chiral symmetry 
breaking. At high enough temperatures, on the other hand, the preferred configuration is 
that of disjoint D8- and D8-branes; hence, chiral symmetry is restored. 

A disadvantage of the Sakai-Sugimoto models is that, despite being models of sponta- 
neous xSB, the order parameter for such a breaking is absent. In these models there is no 
mode in the bulk geometry from which one can extract holographically the order parameter 
for chiral symmetry breaking. Also, related to this issue, one cannot write an explicit mass 
term for the localized fermions because there is no direction transverse to both the color and 
the flavor branes along which to stretch an open string. So far, there are three proposals 
for how to modify the models in order to be able to compute the xSB order parameter and 
incorporate a bare fermion mass: 

1) Open string tachyon: Ref. [3] argued that one should include the dynamics of an open 

string stretched between the flavor branes, whose scalar mode (tachyon) transforms 
as the bifundamental of U(N{) x U(N{) and has the right quantum numbers to be 
holographically dual to the fermion mass and condensate. Refs [12, 13, 14] incorporated 
this open string mode into the system using the tachyon-DBI action [15] and showed 
that the tachyon has a normalizable and a non-normalizable mode, identified with the 
order parameter for %SB and the fermion mass, respectively. 

2) Open Wilson line: In holographic QCD and NJL, one can make a gauge-invariant oper- 

ator out of the left-handed and right-handed fermions, ipiL and ip J R , which are localized 
at different points in the x 4 -direction by including an open Wilson line between them. 
The vev of this open Wilson line (OWL) operator is an order parameter for %SB, ar- 
gued in [1] to be holographically dual to a Euclidean string worldsheet bounded by the 
U-shaped flavor brane configuration. 

3) Technicolor D4-branes : Ref. [16] added a technicolor sector to the original set-up by 

introducing "technicolor" D4-branes parallel to the color D4-branes. The flavor branes 
form a U-shaped configuration in the near-horizon geometry of the technicolor D4- 
branes, and hence chiral symmetry is spontaneously broken in this sector. A string 
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stretched between the technicolor and the color D4-branes is a gauge boson which can 
then mediate the chiral symmetry breaking of the technicolor sector to the QCD (or 
NJL) sector. As a result, mass can be generated for the fermions of the QCD sector 
through a four-Fermi interaction. To calculate the fermion mass, one has to calculate 
the area of a Euclidean worldsheet bounded by the D4, D8 and D8-branes. 

In this paper, we focus on the second of these proposals and use it to compute the xSB 
order parameter as a function of temperature and constant background electric and magnetic 
field. In section 2, we briefly review holographic NJL, along with Aharony and Kutasov's 
proposal for computing the %SB order parameter (OWL). In section 3, we investigate the 
effect of finite temperature on (OWL). We find the surprising result that (OWL) increases 
with temperature as one approaches the chiral phase transition temperature from below, 
even though the dynamically generated fermion mass decreases. Both vanish above the 
transition temperature, as is expected in the chiral symmetry restoration (%SR) phase. We 
explore the possibility that this unexpected behavior of (OWL) in the xSB phase is merely 
the reflection of the open Wilson line inside the operator rather than the fermion insertions 
at the endpoints but find evidence that this is not the case. In fact, we find that there are 
other models, such as the holographic Gross-Neveu and compact Sakai-Sugimoto models, 
for which (OWL) decreases with temperature or is independent of temperature, respectively. 
This suggests that the behavior of (OWL) also strongly depends on the background in which 
the OWL operator is being studied. 

In section 4, we find that (OWL) decreases (increases) with a constant background 
electric (magnetic) field, regardless of whether or not the temperature is turned on. At 
high enough electric field (in some appropriate units), we show that (OWL) vanishes. Such 
behavior is also seen for the usual chiral order parameter as a function of background electric 
or magnetic field in the standard NJL model, and is consistent with the inhibition of the chiral 
transition by a background electric field, as well as the catalysis of %SB by a background 
magnetic field. 

In section 5, we demonstrate the existence of a family of generalized OWLs which cor- 
respond to the string worldsheet curving in additional directions. These different xSB order 
parameters in holographic NJL may correspond to open Wilson line operators with addi- 
tional insertions in the five-dimensional gauge theory description. By analyzing the global 
structure of their holographic dual worldsheets, we show that, as expected, the vevs of these 
chir ally- charged generalized OWL operators vanish in the xSR phase of the model. In sec- 
tion 6, we consider more general OWL operators for which the fermions are located at two 
different points in spacetime, and show that their vevs also vanish in the xSR phase. 

2 Review of OWLs in holographic NJL 

In this section, we review the computation of the x^B order parameter in holographic NJL 
as proposed in [1]. We start the section by first presenting a short review of the holographic 
NJL model [3]. For simplicity, we set N { — 1 in our discussion as the generalization to 
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N { <C N c is straightforward. 



2.1 Holographic NJL 

Consider the non-compact Sakai-Sugimoto brane setup described in the introduction, where 
the flavor D8- and D8-branes are separated by a coordinate distance £q in the non-compact 
x A direction along which the N c D4-branes are extended. In this set-up, the left and right 
fermion modes, ipi and i/jr, are localized at D8-D4 and D8-D4 intersections and so interact via 
D4-brane gauge fields (and scalars) whose strength is controlled by an effective dimensionless 
coupling A cff = X 5 /£ . 

At weak coupling, A c ff <C 1, and in the regime where g s N c is small so that stringy effects 
can be neglected, integrating out the D4-brane gauge fields (and ignoring the interaction 
due to the exchange of D4-brane scalars) generates an effective action for the fermions with 
a non-local interaction term [3] 

Sint~#5 J d 4 xd 4 yG(x -y,£ ) ip[{x) ■ tp R (y) ^ R {y) ■ ip L (x) , (1) 

where G(x — y, £$) is a 4 + 1-dimensional scalar propagator and the dot in the parentheses 
denotes contraction over color indices. Thus, the theory of the fermions is a non-local version 
of the NJL model with a natural cut-off of A = i~ l . At weak coupling, the analysis of the 
gap equation shows [3] that ip\ j (x)ipFi(y) develops a condensate breaking the chiral symmetry 
even at arbitrarily small values of \ c g. This is different from the standard NJL model where 
the condensate cannot form below a critical coupling. 

For A e fj ~ 1, the single gluon exchange approximation that leads to (1) is not valid. At 
strong coupling (A c fj ^> 1), a better description of the system is by gauge-gravity duality. 
In the probe approximation, one considers the flavor branes in the near-horizon geometry 
created by N c extremal D4-branes given by the metric 

ds 2 = (|) 3/2 (dt 2 + + Q ~ 3/2 {du 2 + u 2 dti£j , (2) 

where d£l\ is the metric of a unit 4-sphere. The characteristic length parameter R is given 
by 

R 3 = ng s N c l 3 s = ^g 2 N c l 2 s: (3) 
where g s is the string coupling. Also, the dilaton and the 4-form RR-flux F( 4 ) are given by 

a /«\ 3/4 r, 2ttA" c 

e * = 9s{x) , ^(4) = -^e (4) , (4) 

with V4 and e( 4 ) being the volume and the volume form of the unit 4-sphere, respectively. 
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The dynamics of a D8-brane (and a D8-brane) is determined by its DBI plus Chern- 
Simons action. Solving for the equations of motion for the gauge field, one can safely set the 
gauge field equal to zero and just work with the DBI part of the action 

Sdbi = -fJ-s J d 9 a e^^-det(g ab ), (5) 

where /x 8 is a constant and g ab = GMNd a x M d b x N is the induced metric on the worldvolume 
of the D8-brane, with Gmn being the metric for the near horizon geometry of the color 
D4-branes. Since we would like to determine the shape of the flavor branes as a function 
of the radial coordinate u, we choose an embedding where the flavor branes form a curve 
u = u(x 4 ) in the (u, x 4 )-plane, extend in M 3 ' 1 x S* 4 and are subject to the boundary condition 

w(±y)=w A ^oo. (6) 

The DBI action (5) then reads 

S Dm = -C J d^x dx A «yi+(^)"V, (7) 

where C = ^V^/g s and the primes here denote the derivative with respect to x A . 
The first integral of the equation of motion for u(x 4 ) is 

= < (8) 



where uo parametrizes the solutions. The simplest solution compatible with the boundary 
condition (6) is x A = ±£q/2 which represents disjoint D8- and D8-branes. Note that this 
"trivial" solution is obtained by setting uq = in (8). For u ^ 0, solving for u yields 

— i(l)> 8 -«4 <»> 

which shows that u has a turning point at uq. This class of curved flavor brane solutions 
labeled by Uq are interpreted as representing the breaking of chiral symmetry: at large u 
the U(Nf) x U(N{) symmetry is manifest whereas at u — u , where the branes join, there is 
just one U(N{) factor. For a fixed £ , the curved solution is more energetically favored over 
the disjoint solution indicating that holographic NJL at strong coupling has a vacuum with 
XSB. 

Integrating (9) results in 

r e °/ 2 C>3/2 , q I . 

/ & = (10 ) 

where we defined the convenient set of integrals 

T [°° du 1 _i( n+6 ) /n + 6 1\ . . 

I n — I — ; — = -u n B[ ,- , n > 0. (11) 

Juo WV« 8 - «o 8 8 V 16 '2/' 1 ; 
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Equation (10) implies that for each £ there exists a unique solution representing holograph- 
ically a %SB phase of the dual gauge theory. It also indicates that, the larger the asymptotic 
separation, the lower the x^B scale Uq. Although equation (10) by itself does not put bounds 
on £ , other considerations do so. As u — > 0, one enters a regime of high curvature for which 
the supergravity approximation is no longer valid, putting an upper bound on £ . Also, in 
order to ensure the stability (or metastability) of the flavor branes at large u, one has to 
require £ ^> £ s which gives a lower bound on £ . 



2.2 OWLs as order parameters for xSB 

Since ipi and ipn are localized at different points on x 4 , a Wilson line must be inserted in 
the standard (V'l xSB order parameter in order to render it gauge invariant. Thus, we 
define the OWL operator 



r r £o/2 i 

OWL^x") = ^(x»,x 4 = -£ /2)V exp / {iA A + <$>)dx 4 ^ iR (x^,x 4 = £ /2), (12) 



where A 4 is the component of the gauge field in the x 4 -direction, and $ is a scalar. Note 
that when ipn and ip R are weakly-coupled to and $, the operator defined in (12) reduces 
to the usual order parameter for chiral symmetry breaking in the NJL model. 

One can make many more non-local operators of the above type by choosing different 
contours for the Wilson line, or by inserting other operators along the contour. These 
operators are interesting in their own right, and we will discuss some of them in later sections. 
For now, we focus on the simplest such operator, given by (12). 

Analogous to holographic Wilson loops, it was proposed in [1] that at strong coupling 
the operator (12) is holographically dual to a Euclidean worldsheet bounded by the contour 
on which the operator is defined at u = u\ and the flavor branes, and localized in the rest 
of the directions. Furthermore, the one-point function of OWL^i^) is given by 

(OWL]) ~ 5{e~ SF , (13) 

where Sp is the regularized action of a Euclidean fundamental string whose worldsheet was 
described above. The leading contribution to the string action is given by the Nambu-Goto 
action 

S F = / da 1 da 2 J^etha^, with h a p = G M Nd a x M dpx N . (14) 

Zna J 

One can check that the worldsheet described above is a solution to the equations of motion 
arising from the Nambu-Goto action. 

To calculate the regularized area (action) of this worldsheet we first choose the gauge 
x 4 = a 1 and u = a 2 , and substitute this back into (14) to obtain 

Sf = 2~W J d A^-u(x% (15) 
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where u\ is the cutoff. By changing the measure of the integral to du, we obtain 




For a finite result, the linear divergent piece in (16) can either be absorbed within (OWL]) 
or subtracted away by a Legendre transform of the action (16) with respect to u [17]. Using 
(3) and (10) together with A 5 = (2ir) 2 g s N c £ s: one arrives at 

S F = -cA cff , (17) 

where 

c=—^b(— ,-)b(—,-) « 0.008. (18) 

128tt 2 \16 2/ V16 2/ v ; 

Thus, (13) yields 

(OWL]) ~ 8je cX ° s . (19) 

The above result represents the leading contribution in A e g to (OWL]). The next-to-leading 
contributions come from two sources: the coupling to the dilaton in the worldsheet action 
[18], and the fluctuation determinant around the saddlepoint. For the worldsheet in our 
discussion, these effects have been calculated in [19]. In this paper, we only consider the 
leading contribution in A cff to (OWL]). 

Note that the holographic computation leading to (19) is valid for A e g 1, which implies 
that (OWL]) is exponentially large. The exponential behavior of (OWL]) is unexpected 
from our experience with QCD, or the NJL model where the order parameter is related to 
the dynamically generated fermion mass. The unexpected behavior of (OWL]) has been 
attributed in [1] to the existence of the open Wilson line in the definition of the operator in 
(12), not the fermion bilinears. In what follows, we will see more surprises in the behavior 
of (OWL]) when we consider the model at different external conditions. 

3 OWLs in thermal holographic NJL 

3.1 Holographic NJL at finite temperature 

To consider the holographic setup at finite temperature, one studies the dynamics of the 
flavor D8- and D8-branes in the near-horizon geometry of a stack of non-extremal D4-branes 
with the metric 

^ = (in - ^ + **■) + (ir /2 (^ + - 2 ^)' /M = * - < 2 °> 
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where ut is the horizon radius, which is related to the inverse temperature (3 of the black 
brane by 



ut = - — . (21) 



16tt^R 3 

9 (3 

The dilaton and the 4-form RR-flux are the same as in the previous section. 



In order to determine the vacuum, we use the same embedding and boundary condition 
as in the previous section. Hence, the DBI action for a D8-brane in this geometry is given 
by 



'DBI 



-C J d 3+1 x dx 4 w 4 ^//(w) + (|) V. (22) 



The equation of motion gives 

,2 1 fU\ 3 

which upon integration gives 



m ' 2= ^Q) 3/(m)[m8/(m) " m ° 8] ' (23) 



^0 = ^3/2^4 f°° du ( 24) 

2 °J Ut M 3/2 v / /(m)K/(m) _ m 8 ] ' 

where the turning point u t is given by the largest real root of the equation 

u 8 t - u%u b t - ul = 0. (25) 

The trivial case of u = represents parallel flavor branes descending all the way down 
to the horizon u = Ut- This is interpreted as a phase of holographic NJL in which chiral 
symmetry is unbroken. A given non-vanishing value of u represents a U-shaped flavor brane 
solution. However, for a given value of £q at a particular temperature the analysis of (24) 
and (25) shows that there are two different values of Ut- Thus, there are two branches of 
solutions for which the flavor branes are smoothly joined at some turning point u t . Figure 
1 (a) shows plots of £o versus u t for different temperatures where, for convenience, we have 
fixed the asymptotic separation of the flavor branes to £ = 1. We have also set it! = 1. We 
will refer to the solution with larger (smaller) turning point u t as the short (long) solution. 
For > 0.167, there exists no U-shaped solution for which £q — 1. For this temperature, 
Figure 1 (b) shows that u t — UT does not vanish, indicating that the U-shaped solution ceases 
to exist before falling into the horizon. 

Both short and long solutions are realizations of chiral symmetry breaking, though the 
constituent dynamical fermion mass (u t — UT)/27ia' is different for the two solutions. Energy 
analysis shows that, below a critical temperature of /3~g B w 0.15 (in units of £^ 1 ), the short 
solution is favored over the solution with parallel branes (and also over the long solution); 
hence, chiral symmetry is broken. Beyond this critical temperature, the parallel solution 
is energetically favored so the vacuum represents a chiral symmetry-restored (xSR) phase. 
Note that from (25) it is easy to see that u t > u for non-vanishing w T . For fixed £ , the 
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turning point increases as the horizon radius increases, but does so at a slower rate such that 
the two radii coalesce at high temperature. Figure 1 (b) shows that u t — ut decreases as 
temperature increases (up to the phase transition temperature £C SB ~ 0.15). Although the 
plot is for a specific value of £q — 1, its qualitative behavior stays the same for all £o's. In QCD 
or in the NJL model, if the dynamically generated fermion mass decreases as a function of 
temperature, one concludes that the order parameter ipm) should also decrease. However, 
as we mentioned in the previous section, for the model under discussion there is no particular 
relationship between the dynamically generated fermion mass (u t — u^jl'Kol and the order 
parameter (OWL^). Thus, it is a matter of calculation to see how (OWL^) behaves as a 
function of temperature. 



(a) (b) 

Uf—Uj, 



2.5 - 0.6 - 




Figure 1: Plot (a) shows £o versus ut for various temperatures: = 0.08 (blue), 0.1 (red), 0.12 (green), 
0.15 (violet) and 0.16 (yellow). For convenience, we set R = 1 and cut the plots by the line £q = 1 in order 
to better show the values of u t for which £o — 1. For /3 _1 > 0.167, there exists no U-shaped flavor branes 
with an asymptotic separation of £q = 1. Note that the behavior shown in the plots stays qualitatively the 
same for all £o's. Plot (b) shows that Ut — «t decreases with temperature. The plot is for £q = 1 but the 
qualitative behavior is the same for other values of £q as well. As can be seen, the plot stops at /3" 1 « 0.167 
(beyond which there are no U-shaped solutions) but before reaching this temperature the system undergoes 
a xSP phase transition at P~g S ~ 0.15. 



3.2 OWLs in the xSB phase 

Let (OWL^)^ denote the order parameter at finite inverse temperature (3. In the chirally- 
broken phase of the model where the flavor branes smoothly join at a radial point above 
the horizon, one can show that there always exists a Euclidean worldsheet extended in the 
(u, x 4 )-plane bounded by the flavor branes, and at a fixed point in the other directions. In 
order to calculate (OWL^)^, one has to find the regulated area of the Euclidean worldsheet in 
the black hole geometry (20). We choose the same worldsheet embedding as before: x 4 = a 1 
and u = a 2 . In this gauge, the action of the Euclidean worldsheet becomes 

= -K [ dx* I" (26) 

2™ J Ju( x 4) \/f(u) 

Although the above expression cannot be integrated in closed form, we can provide an 
analytic expression at low temperatures (compared to ^q 1 ), u t ^> Mr- From (25) one can 
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easily see that u t ^> ut implies that u ^> u T , where u is the minimum radius of the U- 
shaped flavor branes at zero temperature. Comparing (26) with (15) shows that (OWL^)^ 
should take the form 



(OWLi) p = 5jexp\c\ eS J2 



m=0 



(27) 



where a m > are to be determined by evaluating (26). Note that a — 1. The above 
expression can be further simplified. Since f(u) = 1—u^/u 3 , the integral (26) indicates that 
only a§k in (27) are non- vanishing where k G Z + . Therefore (OWL^)^ takes the general form 



k=0 



(28) 



In the following expressions, we will include the leading finite temperature corrections to the 
worldsheet action, namely we will keep terms up to the order Uj./uq. Thus, (26) becomes 



2iiat 
to 



dx 4 



"A 



2ira 



jUa 1 - 



u(x 4 ) 
Uj. 



) - ^ ,2 < 



Uf 



1 + 



du 



4^ 3 J M l/2 V / M 8_ M 3 M 5_ M 8' 



(29) 



where 



u t ~ u ( 1 + 



My 

8^ 



The following change of variable 



u — v\ 1 + 



enables one to write the integral expression in (29) as 

to ' "* 



(30) 



(31) 



(32) 



Digression on the Legendre-transformed string action. One may wonder whether a 
Legendre transform of the action along the lines of [17] will take care of the dependence of 
the worldsheet action on u\. Although performing the Legendre transform with respect to 
u gets rid of the linear divergent term in (32) [20], it cannot cancel the subleading terms in 
u\, and in fact introduces new counter-terms to the action. Define 

Sp — > Sf — J dx A yJ det(/i Q( g) tc u u, (33) 

where the integral is over the boundary of the worldsheet, h a p is the boundary metric and 

n u = 7 — F G uu n^d lx u, (34) 
27TQ; 
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is the momentum pulled back to the boundary and n M is a normal outgoing vector to the 
boundary. Calculating det(/i Q/ g) and n u at u = u\, we obtain 

det(M = (^J , ^ = w/mIiJ • (35) 

Substituting (35) into (33) and expanding f(u\) for small temperatures, we arrive at 

showing that the Legendre transform of the action removes the linear divergent piece (but 
leaves dependence on the cutoff in non-divergent terms). 



Nevertheless, in the «a — > oo limit, one arrives at an unambiguous result for the action 



S F = -^R 3/2 4(h-^I 7 ). 
na V 16 / 



At low temperatures, £q can also be approximated as 



to 
2 



R^4(l 3 + ^I 9 ). 



(37) 



(38) 



Expressing u in terms of £ and substituting the result in (37), one obtains 



SV = — cA 



cff 



1 + a T 



where c has already been defined in (18) and 



a,T 



8 7T 

3fi 



Thus, we obtain 



(JL i 

^16' 2 



(OWL]) 



5 Tfi> 



9_ 1 

16' 2 



TDI 7_ 1 

1 16' 2 



1 + 1.09A cff 



4\ 6 



136.2. 



(39) 



(40) 



(41) 



In order to demonstrate that this behavior persists as one increases the temperature, we 
numerically solve for (OWL^)^. For convenience, we set R = 1, i = 1 and 2na' = 1, resulting 
in A 5 = 8/T 2 . In order for the IR sensitivity of the OWL not to be smeared out by integrating 
too far in the UV region, we choose the cutoff ua = 10. As we vary the temperature, we 
use a shooting algorithm to find the values of u t which correspond to keeping £ — 1. Figure 
2 not only shows that (OWLj)^ increases with temperature but also that the transition to 
the chirally-symmetric phase of the model is first order. One can check that behavior of 
the plot stays qualitatively the same for other values of £ as well. We have cut the plot at 
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{OWL) f /{OWL) 




0.98 1 ' ' ' B" 1 

0.05 0.1 0.15 

Fi gure 2: (OWL^)^/(OWL^) increases with temperature in the chirally- broken phase of the model, and 
the chiral restoration transition is first order. The plot is for £q = 1 and, for convenience, we set R = 1, 
2-kq,' = 1 and ua = 10. The plot stays qualitatively the same for larger values of %. 

/3~g B ~ 0.15 where the chiral restoration phase transition happens [11, 21] but, in principle, 
one can continue the plot beyond /3~g B up to < 0.167 which is the critical temperature 
(for a fixed £q = 1) beyond which no U-shaped solutions exist. 

Our analytic and numerical results in this section both show that, for holographic NJL 
at strong coupling and in the the phase, (OWLj)^ increases monotonically with temper- 
ature. This behavior is quite different from what one observes in QCD and in the standard 
NJL model. Indeed, in the standard NJL model at finite temperature, the analysis of the 
gap equation for the order parameter (^iipm) shows that it decreases with temperature. 
Plotting the effective potential for the order parameter Kff((V'x Vir}) versus (ip^ipm) for 
various temperatures, one observes that the absolute minimum of the plot moves towards 
smaller values of (ip^ipm) as the temperature is increased (see for example Figure 1 in [22]). 
At a critical temperature, the minimum is at (^l V'ift) — 0, for which the model goes to the 
chirally-restored phase. 

3.3 Contour-cancelled x^B order parameters 

Although (OWL^) goes over to the usual order parameter (ip^ipiR) at weak coupling for the 
holographic NJL model, its unusual temperature dependence raises the question of whether 
there exist other gauge-invariant chirally-charged operators similar to (12) whose vev shows 
the behavior expected of an order parameter in the NJL model. In fact, there are infinitely 
many OWL operators which can act as chiral symmetry-breaking order parameters. For 
instance, the shape of the Wilson line contour can be changed and other local operators can 
be inserted at points along its length. (Some examples of these modified order parameters 
will be examined in detail in sections 5 and 6, below.) 

Here, let us try to use this freedom to define other order parameters to test one possible 
explanation (already mentioned in [1]) for the unexpected behavior of (OWL?) both at zero 
and at finite temperature: namely, that the value of its vev is dominated by the open Wilson 
line contour inside the operator rather than the fermion insertions at the end points. A 
way of testing this is to consider ratios of operators with different Wilson line contours. In 
particular, an appropriate ratio could (approximately) cancel the contour dependence of the 
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OWL, leading to a %SB order parameter whose vev is dominated by the fermion insertions. 
So, does such a contour-cancelled order parameter decrease with increasing temperature? We 
will now argue that it does not. This suggests that the unusual behavior of the OWL vev is 
not due to its being dominated by the Wilson line as opposed to the fermion insertions, and 
that the temperature-dependence of the OWL order parameter reflects a physical property 
of the xSB transition in holographic NJL, and is not just an artifact of a poor choice of order 
parameter. 

One way to cancel the contour dependence is to divide the OWL operator by the square 
root of the vev of the closed Wilson loop (CWL) along a contour that traverses the open 
contour in one direction then returns along the same contour to its starting point. However 
this clearly does not work because the CWL vev has no temperature dependence as the area 
enclosed by the contour vanishes. A better choice is to open up the CWL contour to reduce 
such cancellations. For example, a particularly simple and symmetrical choice is to choose 
the CWL to be a circle in the x 3 -x A plane of radius £q/2, and to choose the OWL contour 
to be the semicircle of the same radius between the fermion insertions. In computing the 
saddlepoint contribution, we consider the corresponding worldsheet for the CWL to have 
boundary on the circle at u — U\, while the worldsheet for the OWL has boundaries on 
the semicircle at u = and on the D8-branes. These worldsheets are roughly depicted in 
Figure 3. Then we can construct our contour-cancelled order parameter as 



u OWL CWL 




Figure 3: The shaded areas represent string worldsheets ending on a semi-circular open Wilson line in the 
xz-xa plane u = u\ (green) and on a circular Wilson loop of the same radius in the same plane (blue). For 
the open Wilson line, the X's correspond to the endpoint sources and the worldsheet's dashed boundary lies 
along the D8-brane, which is not shown. u± is the altitude of the D8 brane turning point, while u s is the 
altitude of the turning point of the circular Wilson loop's string worldsheet. 



OWL J 

;= uvv ^ (42) 
V(CWL) 

so that 

(Of) ~ ^e~ (SF °p-"^ 5F ciosod), (43) 

where Sp is the action of the open string worldsheet with the semi-circular contour and 
D8-branes as its boundary, and Sp zloscd is the action of the worldsheet with the circular 
Wilson loop as its boundary. 
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Note that this kind of contour-cancelled order parameter has the advantage of automat- 
ically being UV regulated, since the contributions from the string worldsheets near u — u\ 
cancel. A disadvantage of this type of order parameter, however, is that having the contour 
bend in some spacetime direction (a; 3 in the figure) necessarily breaks space-time rotational 
invariance, which was preserved by Aharony and Kutasov's OWL order parameter proposal. 

It is difficult to analytically evaluate the dependence of this contour-cancelled order pa- 
rameter on the temperature. For example, the equation governing the string worldsheet with 
boundary on the CWL is derived from the Nambu-Goto action. Choosing polar coordinates 
(r, 0) in the x 3 -x 4 plane and assuming that the worldsheet is ^-independent so that u = u(r), 
the Nambu-Goto action is proportional to 



where we have rescaled r — > R 3 ^ 2 r to remove extraneous factors, so r = (£o/2)R 3 I 2 . This 
gives the equation of motion 

= 2ru(u 3 - u 3 T )u" + 2uu' 3 - 3r(2w 3 - u 3 T )u' 2 + 2u{u 3 - u 3 T )u' - 3r(u 3 - u 3 T ) 2 , 

with boundary conditions u' = at r = and u = «a at r = r . This has to be solved 
numerically. 

The general result is that the CWL string worldsheet has a much smaller area than the 
OWL worldsheet, and so that ^(CWL) > (OWL). This is reflected in the fact that the 
altitude u s of the CWL string worldsheet is much higher than the turning point u t of the 
D8-brane. This is simple to understand qualitatively: the fact that the D8-brane is extended 
in more dimensions than the string means that its area scales as a higher power of u, and 
so decreases its area by descending to lower values of u. Thus, the the turning point u t of 
the D8-brane is much more sensitive to changes of the temperature (and, therefore to the 
location of the horizon ut) than is the turning point u s of the CWL string. For example, for 
£o = 0.5 (in units of R), u t varies approximately between 2.10 and 2.25 as the temperature 
increases from to the x^R transition, while u s varies only from 38.772 to 38.774 over the 
same temperatures, which is only about 1/75 of the first variation. 

Since the area of the worldsheets is essentially proportional to £o(u\ — u s ) or £o( m a — u t) 
for the CWL or OWL worldsheets, respectively, it follows that the temperature dependence 
of u t completely dominates that of u s . Thus, the exponential increase of the order parameter 
with temperature is hardly any different between the original (OWL) order parameter, and 
the contour-corrected one, (O). 

3.4 xSB order parameter in holographic Gross-Neveu 

Although (OWL]) does not show the expected behavior of an order parameter in the %SB 
phase of holographic NJL, it may well be the case that it shows the expected behavior (that is 
to say, it decreases with temperature) for other intersecting brane models of chiral symmetry 
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breaking and restoration. To put it another way, (OWL^) may not only strongly depend on 
the choice of a contour, but also on the background in which it is embedded. 

For instance, consider a D2-D8-D8 system where N { flavor D8- and D8-branes intersect 
N c color D2-branes at two (l+l)-dimensional intersections, and are separated by a distance 
£q, say, in the y-direction. The massless degrees of freedom at the intersections are left- 
handed and right-handed Weyl fermions which interact via the exchange of gauge fields (and 
scalars) of the color theory whose strength is given by an effective coupling A cff = \ 3 £ . At 
small effective coupling, after integrating out the gauge fields, one arrives at a generalized 
Gross- Neveu (GN) model with a non-local four- Fermi interaction [23] 1 . At strong coupling, 
one studies the dynamics of the flavor branes in the near-horizon geometry of extremal D2- 
branes. This theory, which we call holographic GN, like holographic NJL, breaks chiral 
symmetry. It also shows a x^R transition at a finite temperature. 

To study the model at finite temperature and large coupling and in the probe approxima- 
tion, one considers the flavor branes in the near-horizon geometry of non-extremal D2-branes, 
which is described by the metric 



ds 2 

' '13 



where 



it* = 27^, h(u) = 1 - ^. (45) 



The temperature of the black brane is related to ut by 

1 5 /w r \ 5 / 2 1 

The dilaton and the 4-form RR-flux of the solution take the forms 



(46) 



, u \ " 5/4 2ttN c , x 



where V 6 and e( 6 ) are the volume and the volume form of the unit 6-sphere, respectively. 

For a D8-brane which wraps the 6-sphere and forms a curve u = u{y), the DBI action 
reads 

S DBl = -^Vol(5 6 )^ /2 J d l+1 x dy u'/^h^+^Y'u' 2 . (48) 
The equation of motion gives 

= u'\ (49) 



1 There is another holographic realization of the GN model based on intersecting D4-D6-D6-branes which 
gives rise to a non-local four- Fermi interaction [24]. 
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where, as usual, u = is a parallel brane solution which is not energetically favored at low 
temperatures. For u ^ 0, we have (at small temperatures) 



where 



to 
2 



du 



, w 5 / 2 ^h(u)[u 7 — u 2 u\ — Mq] ' 



i? 5 / 2 M f(^5-^15), 



1 -i(n+5)„/n + 5 1 



Uo WVu 7 - mo 7 7 



'5 



/n + l\ 



n > 0. 



(50) 



(51) 



The action of a Euclidean worldsheet bounded by the flavor branes, stretched in the 
(u, y)-plane and point-like in other directions is given by 



Sv = 



1 



2na 

to 



dx 



U A 



du 



u(x 4 ) \Jh(u) ' 



K 2^7- ('-W-^'^-f^' (52) 

where in the second line we have only kept the leading term in the small temperature 
approximation. Eliminating m between (50) and (52), dropping the -UA-dependent terms in 
(52), and keeping the leading temperature-dependent term, we get 



S F = -6(Acfr) 1/3 
where b and bx are given by 



1 h ( £0 Y 0/3 



(53) 



b = I 



1 1/3 



28 



7 H?.i)J 



/4 1\ 
B( f>5 )»0.63, 



br = — 



1 

5~6 



-i 10/3 r 



14 7T 



/S(f,l)_ 



23- 



B 



10 1 

7 ' 2 



5 



4 1 
7' 2 



5 



5 1 
7' 2 



10.5. 



(54) 



So, we have 



" " cff ' 



(55) 



showing that the (OWL^)^ decreases with temperature in the small temperature approxi- 
mation. 

This behavior persists for higher temperatures, as shown by numerical calculations. For 
example, in Figure 4 we show the results of such calculations where we have set R% — 1, 
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2na' = 1 and u\ = 10, and plotted (OWL])^/ (OWL]), as well as u t — ux, versus temperature 
for U-shaped branes with an asymptotic separation of £q = 1. With our choice of 
parameters, A e g = 47r, indicating that holographic GN is in the strong-coupling regime. 
Note that, although in the plots of Figure 4 we considered temperatures up to ~ 0.22 
beyond which there is no U-shaped solution with £ = 1, the system undergoes a xSR phase 
transition at /3~ SB ~ 0.2. Figure 4 (b) shows that for temperatures up to P^SB' (OWL*)/? 
decreases with temperature. Also, since (OWL])^ does not vanish at /3~g B , the xSR phase 
transition is first order. The qualitative behavior of the plots presented in Figure 4 is 
independent of the particular choice of the cut-off u\ we have made, and also stays the same 
for all allowed values of £q. 



(a) (b) 
u t -u T (OWL)J(OWL) 




0.05 0.1 0.15 0.2 0.05 0.1 0.15 0.2 



Figure 4: For U-shaped branes for which £o = 1, we set R3 — 1, 2ira' ~ 1 and ua — 10, and numerically 
plotted (a) u t — ut, and (b) (OWL]) ,3 /(OWL]) versus temperature The plots show that, in the xSB 

phase of holographic GN, both ut — ut and (OWL^)^ decrease with temperature , and the transition is first 
order. 



3.5 xSB order parameter in the (compact) Sakai-Sugimoto model 

As another example, consider the Sakai-Sugimoto model, where x 4 is now a circle with radius 
Ro and the flavor branes are located at the antipodal points of the circle, £0 = ttRq. At finite 
temperature, there are two competing background geometries from which one determines the 
deconfinement temperature The deconfinement temperature is obtained when the two 
geometries are equally energetically favorable. At low temperatures when the model is in 
the confined phase, the background geometry is the same as the zero temperature geometry 
[2] but with Euclidean time periodically identified. Also, in this phase, the flavor D8- and 
D8-branes smoothly join at the tip of the background geometry, thereby realizing chiral 
symmetry breaking. At high temperatures, the model is in the deconfined phase and the 
background is identical to (20) except that the x 4 -direction is a circle. In this phase, the 
preferred configuration is that of disjoint flavor branes. Hence, chiral symmetry is restored. 
One can easily show that, in the confined phase where chiral symmetry is also broken, 
(OWLj)^ does not change with temperature and equals its value at zero temperature [1] 

(OWLf)^ = (OWL]') ~ 8ie Xs/187rRo , V /T 1 < (56) 
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This is the expected behavior of the %SB order parameter (ip^tpiR) in the confined phase of 
finite temperature QCD [25]. Also, (OWL^)^ vanishes in the deconfined phase (see section 
3.6 for more details). 

For the generalized Sakai-Sugimoto model where £q < nR, there also exists an interme- 
diate phase [10] for which the theory is deconfined while chiral symmetry is broken. For this 
phase, an analysis similar to what we did in this section shows that (OWL^)^ increases with 
temperature up to the xSR temperature. 



3.6 OWLs in the xSR phase 

In the high-temperature phase (temperatures large compared to i^ 1 ) of the model, the flavor 
branes are parallel. There exists a Euclidean worldsheet, bounded by the parallel D8- and 
D8-branes, which penetrates into the horizon. In order to better understand the behavior 
of this worldsheet, especially inside the horizon, we use Kruskal coordinates. This will be 
useful later in our discussion of modified x^B order parameters in section 5. 

To reduce clutter, we scale out the coordinates of the metric by replacing 

u — > ut m, x" -> R 3/2 u T 1/2 x", (57) 
so that the background metric (20) takes the form 
ds 2 



IT = « 3/2 

7 2 



- f(u)dt 2 + df 2 ] + u~ 3/2 [/(u)" 1 ^ 2 + u 2 dn 2 ^ , (58) 



where f(u) = 1 — u 3 , and 7 2 = If-R 3 . Note that in (58), t,x and u are all dimensionless. 



3/3 

First, define the tortoise coordinate r by 

dr w 3 / 2 



du u 3 



(59) 



We choose phases in the solution to (59) so that near u = 1 we have e 3r ~ u — 1. Define the 
infalling and outgoing Kruskal coordinates by 

v ;= +e +3(t+r)/2 ) w ._ _ e -3(t-r)/ 2) (6Q) 

so that the metric becomes 

ds 2 4 u 3 - 1 



7 2 9 u 3 / 2 



-e- 3r dvdw + u 3/2 dx 2 + u l/2 dn\. (61) 



Note that vw = —e 3r , so r and u are functions only of the combination vw. Figure 5 shows 
v w and G vw as functions of u. Note that the singularity (u = 0) is at vw = 1 and the horizon 
(u = 1) is at vw = 0. 

Figure 6 shows the black hole geometry in Kruskal coordinates. The red lines are the 
singularities, the blue dashed lines are the UV cutoff(s), and the v and w axes are the 
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Figure 5: vw(u) (upper, blue curve) and G vw (u) (lower, red curve). 

horizons. In Kruskal coordinates the extended black hole geometry is apparent, including 
the second asymptotic region and the past "white hole" (w < 0). 

The straight string worldsheet descending from the OWL between two straight D8-branes 
is extended along u and x A at fixed values of the other (Schwarzschild-like) coordinates. In 
particular, it is at a fixed value of t, say t = t . This corresponds to the line w = —e~ 3t °v 
in Kruskal coordinates, and is shown as the green line in Figure 6. This line passes into 
the second asymptotic region where it has nowhere to end, as there is no OWL insertion at 
u = u\ there for it to end on. Therefore, in the xSR phase of the model there exists no 
finite-area worldsheet [19] to be holographically dual to OWL;(i' i ); hence, (OWL^) = 0. 




Figure 6: The near-horizon geometry of black D4-branes in Kruskal coordinates. The blue dashed lines 
represent cutoffs and the red lines are the singularities. 

We can also examine the OWL worldsheet in Euclidean signature, where the Euclidean 
time t E is periodically identified with the inverse temperature: t E — t E + (3- The period (3 
is given by (21), so that the horizon becomes a regular point with local polar coordinates 
approximately given by (p, 0) ~ ( \J (3{u — Urjjft , 2-rctE/(3) near the horizon, u = ut- The 
worldsheet descends from the OWL at constant t E , say t E = 0, until it reaches p = 0. At this 
point, it passes through the origin and emerges in the antipodal direction, t E = [3/2. The 
worldsheet continues to arbitrarily large u in this direction, since there is no OWL insertion 
at t E = (3/2 for it to end on. Thus, it again gives zero saddlepoint contribution to (OWL), 
consistent with being in the xSR- phase. 
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4 Holographic NJL OWLs in background electromag- 
netic fields 



In this section, we determine how (OWL]) responds to turning on a background electric or 
magnetic field in holographic NJL. We turn on a U(l) background electromagnetic field by 
gauging the U(l)y part of the U(Nf)y symmetry of the model. We adopt the convention 
of [26] (see also [27, 28]) where the fermions have charge one under this U(l) field; namely, 
A = N f ~ lr TrA where A e U(l) and A G U(N f ). We only consider the case Nf — 1 here. 



4.1 Background electric field 



To add a background electric field to the holographic NJL model in, for example, the re- 
direction, we turn on the A\ component of the gauge field on the flavor branes and choose 
the ansatz 



2na'A 1 = -Et. 



(62) 



For a D8-brane forming a curve x 4 = x 4 (u) in the background of black D4-branes (20), and 
with the gauge field turned on as in (62), the DBI action reads 



S 



DBI 



-/is J da 9 e^^-det(g ab + 2vra'F ab ), 



= -C / d 3+1 x du u 



where F a b = F^ u d a x^dbX l/ is the induced field strength on the D8-branes, C = fJ-sV^/ g s , and 
E is the dimensionless electric field. There is also a Wess-Zumino term 



9 - Nc 
owz — 



24tt 2 



Ju 5 (A) = ^- 2 Jaafaf, 



(64) 



which vanishes in our case. The equation of motion for the brane embedding x A {u) is 

\ 



u 4 f(u) 



/ / \-3 \ V2 

/ 1 W«] E 2 



1 - 



f(u) \ R 



dx 
du 



(65) 



In (63) the first term under the square root is positive, independent of the embedding. 
However, the second term can be negative, rendering the embedding unphysical. Indeed, it 
is the analog of the critical electric field for flat branes and is easy to understand. Consider 
an open string with end points on the flavor branes which is also assumed to lie along the 
x 1 -direction. The tension of this string is equal to its energy per unit length. Assume that 
at u = ma — > oo the string tension is {2 / Ka')~ 1 . Therefore, its effective tension at a radius 
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u equals (2ira')~ 1 a/ f(u) (u/R) 3 / 2 . Now an electric field in the ^-direction will cause the 
string endpoints to move away from one another. If the string is at u < u cr , the electric 
field will overcome the string tension and will pair create. This is where the semi-classical 
description of the flavor brane in terms of its DBI action breaks down. This critical radius 
u CT is where the string tension cancels the electric field and is given by 

ul r = u% + R 3 E 2 , (66) 

reproducing the radius at which the second factor in (63) vanishes. Below this critical radius, 
no physical embedding of the flavor branes can exist. 

Since (66) indicates that the horizon is always below the critical radius, the parallel em- 
bedding of the flavor branes (the chirally-symmetric phase) cannot exist. The holographic 
interpretation of this result is that turning on even a small electric field, E <C 1, in holo- 
graphic NJL will wash out its otherwise would-be chirally-symmetric phase. Of course this 
interpretation is solely based on the DBI analysis of the flavor branes. As the flavor branes 
approach u cr the electric field will break open strings and destabilize the system. Taking the 
corrections into account, it is plausible that a parallel embedding of the flavor branes can 
survive, resulting in a restoration of chiral symmetry in the model. A U-shaped configura- 
tion, on the other hand, can exist as long as the minimum radius u t of the solution stays 
above u cr . 

One consequence of having this critical radius around is that there is a bound on the 
dynamical constituent mass of the fermions. There is a critical radius for an open string 
stretched from the tip of the flavor branes to the horizon of the geometry. The mass (in 
string units) of this string equals u t — ut- Since the flavor branes cannot descend below u CT , 
the mass of the string (the dynamical constituent mass of the fermion) cannot be less than 
u CT — a?. Equivalently, since u t — ut is related to the asymptotic distance £q between the 
flavor branes, one deduces that there is a maximum £ max beyond which there is no U-shaped 
flavor branes. 

To simplify the equation of motion while still capturing the essentials of turning on a 
background electric field, we first consider the zero temperature limit of the equation, thus 
describing holographic NJL at zero temperature but with a finite electric field. 



Zero temperature 

At zero temperature, f(u) = 1, in which case the critical radius becomes u 3 T = R 3 E 2 . The 
= solution in (65), which represents the parallel branes, is not a physical solution for the 
reasons explained above, at least in the DBI approximation employed throughout this paper. 
Thus, we assume u* ^ 0. Note that, in order to trust the calculations, we assume that both 
u t and u cr stay away from the region of high curvature of the background geometry. Setting 
f(u) = 1 in (65) and solving for dx 4 /du, we get 

/d*V = (lY 3 —-$- -, ui = R 3 E 2 . (67) 

\du) \RJ u 8 -u 5 u 3 ci -u* 
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When there is no electric field, the minimum radius for the U-shaped flavor branes (the 
turning point) is at u t = Turning on a small E field, one expects a small deviation in u t 
from u*. Assuming ^- <C 1, we obtain 



To analyze the behavior of the solution near the critical radius, one has to take a different 
limit. As decreases to arbitrarily small but positive values (to exclude = from the 
limit), the turning point approaches the critical radius. To leading order, we have 



Figure 7 (a) shows the number of U-shaped branes with a fixed asymptotic separation of 
£ = 1. For the plot, we set R — l. There is a maximum electric field E mSLX « 0.158, beyond 
which no U-shaped flavor branes with £ Q = 1 can exist. For E < E max , there always exist two 
solutions (similar behavior is also observed for other values of £o). Using the energy argument, 
one can show that the solution with larger u t is always energetically favored. Figure 7 (b) 
shows the plot of u t of the energetically-favored solution (with £ = 1) versus E. Denoting the 
turning point (of the energetically-favored solution) corresponding to the maximum electric 
field by Wmin, we observe that u m m ~ 0.43. Note that u m i n > u cr (E max ) m 0.29. The fact 
that u min > u CT indicates that the E max of the boundary theory (holographic NJL) is not 
related to the maximum electric field on the flavor branes associated with the breaking of 
open strings on the branes at u = u CT . 

Thus, as one increases the electric field, the energetically favored U-shaped brane goes 
further into the bulk until its turning point reaches u inin . At the same time, the turning 
point of the less energetically favored solution increases with the electric field until it, too, 
reaches -u min . The picture emerging is that the two U-shaped solutions approach each other 
as electric field increases, and coalesce when E = E mSuX . 

Since the Euclidean worldsheet from which (OWL^) is calculated does not directly couple 
to the Ai gauge field on the brane, one deduces from Figure 7 (b) that the (regulated) area of 
the worldsheet increases with the electric field. Hence, (OWL^) decreases quadratically with 
the electric field, as long as the electric field is small. This is exactly the behavior observed 
in the standard NJL model for the dependence of the dynamically generated fermion mass 
on a constant background electric field [29]. 

To calculate the leading-order dependence of (OWL|) on the electric field, we first find 
the area of the Euclidean worldsheet. After a bit of algebra, one arrives at 




1 



3 



(68) 




(69) 




(70) 
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Figure 7: (a) £ a versus u t for different values of the electric field: E = 0.01 (blue), 0.05 (red), 0.1 (green) 
and 0.13 (violet). We set R — 1 and cut the plot by the £q = 1 line, since we are interested in the number of 
solutions with an asymptotic separation of £o = 1 • The plot shows that there are two branches of solutions 
approaching each other as E approaches a maximum -E max ~ 0.158, beyond which there is no solution with 
£o = 1- (b) The turning point of the energetically- favored solution decreases with the electric field, reaching 
Itmin ~ 0.43. 



In general, Sp admits the expansion 



S F 



Next, we solve for from 



2na 



sua 



n=0 



U, 



3n 



(71) 



du 



(72) 



Substituting u* back into (70) and removing the linear divergent term from the action, we 
arrive at 



Sf = — cA 



off 
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(73) 



where 
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(74) 



and c is the same as in (18). Thus, (OWL^) reads 

™^l- 0.029A cff ^V 
(OWI4) \R 



(75) 



where (OWL^) # denotes the order parameter for holographic NJL with a background electric 
field turned on. This shows that, in the %SB phase of the model, (OWL.) E decreases with 
electric field, at least for a small electric field. 
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{OWL) E /(OWL) 




Fi gure 8: At zero temperature, (OWL-)b decreases with the electric field E. To obtain the plot, we set £ = 
1, R = 1, 2-Kct 1 = 1, and ma = 10. When the electric field reaches its maximum value of E max w 0.158, the 
worldsheet reaches the radius u t — w m in — 0.43, thereby minimizing (OWl4)s to (OWL-)^ w 0.9 (OWL^). 



For arbitrary values of the electric field up to ^ ~ 0.158, we have solved the equations 
numerically and a plot of (OWI4).e/(OWI4) versus E is given in Figure 8. To obtain the 
plot, we kept £q = 1, R = 1 and 27ra' = 1. Also, in order to avoid smearing out the IR 
sensitivity of (OWL) by integrating too far in the UV region, we chose the cutoff «a = 10. 
Again, the qualitative behavior of the plot is not sensitive to the particular choice of the 
parameters we have made. The result obtained in Figure 8 qualitatively agrees with field 
theoretic computations in the standard NJL model (see for example [29]), where turning on 
a constant background electric field has been shown to inhibit chiral symmetry breaking. 



Finite temperature 

Having done the analysis for zero temperature, the computation at finite temperature is 
almost the same. We first do the computations when both electric field and temperature are 
small: E < 1 and £ //3 < 1. From (65), 

bH =7^ p 1 5~1 8> < = u t + K v , 76 

\ du / j{u)\RJ v? — u b u 6 cr — ui 

where u cr is given in (67). To ensure the reality of the solution, the = branch is excluded. 
Since the critical radius u cr is always above the horizon, the turning point is the largest real 
root of u 8 — M 5 M^ r — u 8 , = 0. To leading order in temperature and electric field, the turning 
point u t is 

"•^•( i+ S)- < 77 > 



Integrating (76) yields 

^ = R 3 / 2 ui 



du 



u 3 / 2 \J f{u)^u & — u h u\ T — 



u 



8' 



R 3/2 ut (h + ^u 3 T I 9 - ^u 3 CI I 9 ) . (7* 
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Again, one can show that a Euclidean worldsheet extended in the (a; 4 , w)-plane, bounded by 
the flavor branes, and fixed in the other directions is a solution to the equation of motion 
obtained from the Nambu-Goto action. The action of this Euclidean worldsheet now becomes 

s F = i l dx < *■ 



27T«' J J u ( x i) y/f(v,y 

^uji -^)- -^R^ut r (1 + - - t 

KR 3/2 <(h-^u 3 T I 7 -^ull 7 ). (79) 



W * v x 16 J ' 16 

where in the last line above, we subtracted the linear divergent term and ignored the sub- 
leading terms in u\. Eliminating u* between (78) and (79), we get 

S F = -cA cff [l + a T (|) - a s (|) 22 2 J , (80) 

where c, and as are given in (18), (40) and (74), respectively. Thus, one obtains 

(OWL?) ett \RJ : 1 ; 

Since the leading terms in and 22 do not mix, (0WI4')J/(0WI4) retains the same 
qualitative behavior that it had when either or 22 vanished. 

We have also done the numerical calculations for arbitrary values of temperature and 
electric field, for an energetically favored U-shaped brane with a fixed asymptotic separation 
of £q = 1. Figure 9(a) shows Ut — ut versus E for fixed values of f3~ l . This shows that 
the maximum electric field 22 max , beyond which the U-shaped solution (with £ = 1) cannot 
exist, decreases with temperature. Figure 9(b) shows u t — ut versus fi^ 1 for fixed values 
of the electric field, which demonstrates that the maximum temperature beyond which our 
U-shaped solution does not exist decreases with the electric field. 

A three-dimensional plot of u t — ut as a function of and E is given in Figure 10(a), 
which shows that u t — ut monotonically decreases with both and E. (OWL|)^/ (OWL^) 
as a function of f3~ l and E is shown in Figure 10(b). Note that (OWL^)^,/(OWL^) mono- 
tonically decreases with E but, unlike u t — ut, it monotonically increases with f3~ l . In both 
directions, the transition to the chirally-restored phase is first order. 

We would like to emphasize that we wanted to know how (OWL^) responds to turning 
on a background electric field in holographic NJL, which is why we took the ansatz (62). 
We could have generalized this ansatz by taking 27ra'A 1 (t,u) = Et + g(u) where, from the 
asymptotic behavior of g(u), one reads off the current (in the x 1 -direction) in the boundary 
theory. One then finds that, in the x^B phase, the energetically favored solution is the one 
with no current turned on. This is to say that the vacuum is in an insulating phase with 
broken chiral symmetry (see [26, 27, 28] for more details). Thus, turning on a current will 
have no bearing on the computations we performed here for calculating (OWLj) in the xSB 
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Figure 9: (a) ut — ut versus E for various temperatures: /3 _1 = (black), 0.11 (blue), 0.13 (red) and 
0.14 (green), (b) ut — ut versus for different values of the electric field: E = (black), 0.1 (blue), 0.12 
(red) and 0.14 (green). The plots in (a) and (b) are for U-shaped solutions with to = 1 (and R = 1) but the 
qualitative behavior is the same for all £q. 

phase of the model. However, turning on a current, one finds that there exists a special 
value of the electric field (fixed in terms of the current) for which a parallel embedding of the 
flavor branes can be physical [26]. Hence, the high-temperature vacuum is a conductor with 
no chiral symmetry breaking. Since the string worldsheet from which we calculated (OWL]) 
does not couple to A\ and is bounded by parallel branes, the general analysis in section 3.4 
for (OWI4) in the high-temperature phase of the model goes through. That is, the string 
worldsheet penetrates the horizon and extends to the second asymptote of the background 
geometry, in which there is no OWL insertion to end on. Therefore, in this chirally-restored 
conducting phase of the model, (OWL^) also vanishes. 



4.2 Background magnetic field 



To turn on a constant background magnetic field in, for example, the x 3 -direction of the 
boundary theory, we choose the following ansatz for the U(l) gauge field on the D8-brane 



2na'A 2 = Hx 1 , 

where H is dimensionless. The DBI action for the D8-brane reads 
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The first integral of the resulting equation of motion for x 4 (u) is 
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52) 



13) 



(84) 



Since both factors inside of the square root are positive, unlike the electric field, the 
magnetic field does not impose a condition for an embedding to be physical. This implies 
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Figure 10: (a) u t -u T as a function of the electric field E and temperature (b) (OWl^/(OWL^) as 
a function of E and /3 _1 . These plots are for U-shaped solutions with Iq = 1 (and R = 1) but the qualitative 
behavior is the same for all £q. 

that both parallel and U-shaped embeddings exist. The u* = branch in (84) represents the 
parallel embedding where holographic NJL is in the chirally-symmetric phase. Since the A2 
gauge field does not couple to the Euclidean worldsheet (which is only extended in the (x 4 , «)- 
plane), it neither changes the equation of motion nor the boundary conditions. Therefore, 
the worldsheet behaves as if there is no magnetic field. As we showed in the previous section, 
such a worldsheet dips inside the horizon and extends to the second asymptotic region of 
the black hole geometry. This implies that this worldsheet cannot be dual to a one point 
function of the OWL operator. Hence, in the chirally-symmetric phase of holographic NJL 
with a constant background magnetic field one must have (OWL^) = 0. 

We will now analyze how a constant background magnetic field affects (OWL^) in the 
xSB phase. 



Zero temperature 

We turn off the temperature by setting f(u) = 1 in (83) and (84). Solving for dx 4 /du gives 

(SMirW^F^r 4= ^ 2 (85) 

For 7^ 0, the embeddings are U-shaped. Figure 11(a) shows that, unlike for the case 
of a background electric field, there is just one U-shaped solution for a given value of the 
magnetic field. As seen from Figure 11(b), the minimum radius u t of a U-shaped brane 
increases with the magnetic field, which implies that the dynamical constituent mass of the 
fermions increases. Note that Ut asymptotes to a fixed value at large H, which indicates that 
the fermion mass does the same. 
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Fi gure 11: (a) £q versus Ut for different values of the magnetic field: H = 0.05 (blue), 0.1 (red), 0.5 (green), 
f (violet) and 10 (yellow). We set R = 1 and cut the plot at the £q = 1 line, in order to magnify the turning 
point of the solution with an asymptotic separation of £o — 1. For each value of the magnetic field, there 
is only one U-shaped solution, (b) The turning point Ut of the solution monotonically increases with the 
magnetic field H, approaching an asymptotic value at large H. 



To see how turning on a background magnetic field affects (OWL^) at zero temperature, 
we first consider the case where the (dimensionless) magnetic field H is small. To leading 
order, the flavor branes have a turning point at 



u t ~ it* 



8\uJ 



(86) 



Also, the integral of (85) can be expressed as 
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2 



R 3 / 2 u 



du 



17) 



Since A 2 does not couple to the worldsheet, the boundary conditions are not affected. Thus, 
the Euclidean worldsheet we considered in the previous sections remains a solution to the 
equations of motion. The area of the worldsheet in the small magnetic field approximation 
is given by 
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where we have dropped the linear divergent term in the third line. Eliminating between 
(87) and (88) results in 
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(90) 



and c is defined in (18). Denoting the order parameter of the model when there is a back- 
ground magnetic field by (OWL-)^, we have 



(QWL|)g 
(OWLf) 



l + 0.029A eff (|)V 



(91) 



Thus, for a small background magnetic field, the order parameter increases quadratically 
with the magnetic field. Note that the above expression can also be obtained by taking 
E 2 —H 2 in (75). 

This result obtained for holographic NJL at strong coupling agrees with the magnetic 
catalysis of the standard NJL model in the presence of a background magnetic field [30] . In 
fact, the authors of [30] have demonstrated that, not only is the critical coupling of the NJL 
model lowered by a background magnetic field, but the condensate increases quadratically in 
the limit of small magnetic field when the coupling is much larger than the critical coupling 
(see equation (58) in [30]). In Figure 12, we have set R = 1, £ = 1, lita' = 1 and «a = 10 
and plotted (OWL|)#/(OWL^) versus magnetic field H for generic values of the magnetic 
field. 

<OWL„>/<OWL> 



Figure 12: At zero temperature, (OWL^/(OWLj) increases with the magnetic field, where for conve- 
nience we set £q = 1, R = 1, 2na' = 1 and u\ = 10. 



Finite temperature 



We have already given a general argument for why, at high enough temperatures for which 
the theory is in the x^R- phase, (OWL^)^ vanishes for any H . On the other hand, in the 
xSB phase with small temperature and small magnetic field, by taking E 2 — > —H 2 in (81) 
we find that 



1 + 1.09A eff (^) 6 + 0.029A eff (92) 



(OWLj) ~ ' uu \P 
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Figure 13: (a) £ versus u t for H = 0.1 and various temperatures: /3 _1 = 0.15 (blue), 0.17 (red), 0.19 
(green) and 0.2 (violet). We have set R=l. (b) £ versus u t for = 0.14 and H = 0.05 (blue), 1 (red), 2 
(green) and 10 (violet). 

Since the leading terms in and H do not mix, (OWI4)£/(OWl4) retains the same 
qualitative behavior that it had when either (3~ l or E vanished. 

We will now use numerical techniques to consider how the vacuum of the model responds 
to simultaneously turning on arbitrary amounts of temperature and magnetic field. From 
the plots in Figure 13, we see that turning on the temperature results in two U-shaped 
solutions, as in the case of vanishing magnetic field. In particular, Figure 13(a) shows £q as 
a function of the turning point u t for H = 0.1 and various temperatures. This demonstrates 
that there is a maximum temperature beyond which U-shaped solutions with £ = 1 cannot 
exist. A comparison of Figure 13 (a) to Figure 1 (a) shows that this maximum temperature 
is greater than the maximum temperature when there was no magnetic field. Figure 13(b) 
shows £q versus Ut for a fixed temperature (/3 _1 = 0.14) and different values of H. Note that 
the qualitative behavior of all of these plots remains the same for all allowed values of £q. 
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Figure 14: (a) ut — ut versus H for various temperatures: = (black), 0.1 (blue), 0.12 (red) and 
0.15 (green), (b) u t — ut versus for H = (black), 0.5 (blue), 1 (red) and 0.10 (green). The maximum 
temperature beyond which there is no U-shaped branes increases slightly with the magnetic field. 



To explore how the turning point of the energetically-favored solution (with £q = 1 and 
R — 1) depends on the temperature and magnetic field, we plotted u t — ut versus H for 
different temperatures in Figure 14 (a), and u t — ut versus for a few values of magnetic 
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Figure 15: (a) u t - u T as a function of H and . (b) (OWL^)^/(OWL^) as a function of H and (3~ l . 

field in Figure 14 (b). Figure 15 (a) shows a three-dimensional plot of Ut — ut as a function 
of and if. This shows that w t — ut decreases with f3~ , while it increases with H. 

The behavior of the (OWL])^ as function of the temperature (of course, only up to the 
%SR temperature) and magnetic field is shown in Figure 15 (b). (OWl4)J monotonically 
increases with both and H . 

5 Generalized OWLs and curved worldsheets 

According to Aharony and Kutasov's proposal [1], in holographic NJL, (OWL?) ~ 5je~ SF 
where S-p is the action of a Euclidean fundamental string whose worldsheet stretches between 
the regularized boundary u = ma and the D8-branes, whose ends are fixed to lie along the 
OWL contour. The string has Dirichlet boundary conditions along the contour at u — 
u\ and Neumann boundary conditions on the D8-branes for u < u\. This worldsheet is, 
furthermore, assumed to lie at constant values of the other coordinate directions. However, 
there exist other string worldsheets satisfying the OWL contour Dirichlet and D8-brane 
Neumann boundary conditions which do not lie at constant transverse coordinate values. The 
OWL and D8-brane boundaries force the string to be extended in the u- and x 4 -directions, 
but we can also allow it to bend in, say, the t- or x 1,2,3 -directions as well. 

In this section we will show the existence of such curved worldsheets and then discuss 
what is the corresponding holographic OWL operator. We will concentrate on the xSR phase 
with parallel D8- and D8-branes, since it is relatively easy to solve for curved worldsheets in 
this configuration. However, in principle, one should be able to do a similar analysis for the 
worldsheets bounded by U-shaped branes in the %SB phase. In what follows, we perform 
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the computations both in Minkowski and Euclidean signatures of the background geometry 
using Schwarzchild or Kruskal coordinates, as appropriate. 



5.1 Lorentzian signature — Schwarzchild coordinates 

We start with Schwarzchild coordinates and work with a form of the background metric in 
which the coordinates are dimensionless, as we did in (58). 

x°-bending solutions 

We look for possible solutions which bend in the x°- (or t-) direction. We choose the gauge 
and embedding 

x 4 = a\ u = a 2 , x° = t(u). (93) 
Then the Nambu-Goto action is 



-,-% rj-c;-'w (94) 



2%a' J V u 3 — 1 

The lagrangian is invariant under constant shifts in t, giving rise to a conserved charge 
density (energy) carried by the string. So the first integral of the equation of motion is 

dt c u 3/2 

(95) 



du («3_l) v / u 3_ 1 + c 2 

where Co is a constant. 

The right side of (95) has an integrable singularity at u 3 = 1 — Cq (corresponding to 
a turning point) and a non-integrable one at the horizon, u — 1 (corresponding to an 
asymptote). Thus, since solutions must start at u — > 1, they all asymptote to the 
horizon, since the turning point is hidden behind it. The only exception is Co = 0, for which 
constant t corresponds to the straight string solution of Aharony and Kutasov. 

For these solutions, the action becomes 

r«A u 3/2 d , 



>F OC 



/ / 3 TT=2 - (96) 

J ^/u 3 - 1 + eg 



This action is integrable at the horizon, turning point and singularity, and is therefore finite 
unless the string goes off to u = oo. In particular, the solutions which asymptote to the 
horizon have finite area, and thus must continue through the horizon in some appropriate 
in-falling coordinate system, like Eddington-Finkelstein coordinates, which are regular near 
the horizon. As this is just a change of variables from Schwarschild coordinates, the equation 
(95) stays the same; the only difference is that the coordinate change allows us to connect it 
smoothly to a solution with u < 1 with the same value of c . Such u < 1 solutions asymptote 
to the horizon and reach the integrable turning point (u 3 = 1 — Cq). What is less clear in 
these coordinates is what happens to the worldsheet after reaching this point. To understand 
this, it is best to go to the maximally extended Kruskal coordinates, which we will do in the 
next subsection. 
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re-bending solutions 



Now let's look for possible solutions which also bend in, say, the x 2 - (or y-) direction. Choose 
the gauge and embedding 

x 4 = a 1 , u — a 2 , x 2 = y(u), (97) 

so that the action becomes 



*^r*- 3/ v^ +M2 « 98) 

The Lagrangian is invariant under constant shifts in y, giving rise to a conserved charge 
density (momentum in the x 2 -direction) carried by the string. So the first integral of the 
equation of motion gives 

^ = ° 2 (99) 

du V /( M 3_ C 2)( M 3_ 1 )' V > 

where c 2 is a constant. When c 2 = 0, this is again the straight string solution of Aharony 
and Kutasov. 

For these solutions, the action is 

Sf k / 77 3 TT7 3 2V ' (100) 
J y/(u 3 - l)(w 3 - 4) 

which is integrable both at the horizon (u = 1) and the turning point (w 3 = c|), except when 
they coincide. For < j C2 j < 1, the string reaches the horizon at finite y and then turns back 
up. Again, Kruskal coordinates are needed to understand the behavior of the string after 
this turning point, and will be discussed below. For | C2 j = 1, the string asymptotes to the 
horizon. Since it has infinite area, the string does not cross the horizon, even in in-falling 
coordinates. For 1 < |c 2 |, the string reaches a minimum radius of u — |c 2 | 2 / 3 at finite y and 
then turns back up. These turning solutions are non-physical (or, infinite-area) since they 
have no OWL insertion to end on. 



5.2 Lorentzian signature — Kruskal coordinates 

The Kruskal coordinates have been defined in (60) in terms of which the background metric 
was given in (61). Here we use them to analyze the x°- and f-bending worldsheets. However, 
first recall that the straight (Aharony-Kutasov) worldsheet which was the main focus of the 
previous sections corresponds to the line w = — e~ 3t °v . This worldsheet, depicted by the 
green line in Figure 6, passes into the second asymptotic region, where it has nowhere to 
end since there is no OWL insertion at u — u\ there for it to end on. Such a worldsheet 
presumably would contribute to a finite temperature two-point function of OWLs. 
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x°-bending solutions 

Choose the gauge and embedding 

x 4 = cr 1 , vw = a 2 := a, v/w = r(a), (101) 
so that the action is given by 

Sp = — — - [ da\ / 

Here, a\ is related to u\ by — e 3r ^ A ^ := cr A , and is negative. Since u and r are definite 
functions of vw, the quantity in square brackets is a definite (known) function of a . In fact, 
it depends only mildly on a, as shown in Figure 16. 
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Figure 16: The square bracket in the fundamental string Lagrangian as a function of a. 



In Schwarschild coordinates, the Lagrangian is invariant under time translations: t — > 
t + constant. This translates to a scaling invariance under r — » r • e - 3Cons t an t (since 
r = — e~ 3 *). The first integral of the corresponding equation of motion gives 

= , ~ C ° ( 103 ) 
rda avV-l + c 2 ,) 

As a check, upon changing variables from a to w, the above expression becomes identical to 
(95). Thus, we can simply integrate (95) (on both sides of the horizon and of the turning 
point, if the string gets there) and join the solutions by changing to Kruskal coordinates 
using the above formulae. This needs to be done numerically, giving the result shown in 
Figure 17. The fact that the worldsheets look like straight lines is misleading. In actuallity, 
they deviate very slightly from straight lines, essentially by the logarithm of the function 
shown in Figure 16. If one had chosen a larger value of u\ (that is, larger than the choice of 
ma = 1 which is shown), then the curvature would be apparent. There is no good reason for 
ending the |c | < 1 lines at u = u\ in the second asymptotic region. The c = 1 line is just 
tangent to the singularity. The Cq — > oo line becomes lightlike 2 . 



2 It is not clear that worldsheets which hit the singularity (region of high curvature) can be trusted. It may 
be possible for such worldsheets to avoid the singularity by following complexified geodesies to the second 
asymptotic region [31]. 
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Figure 17: The gray lines are string worldsheets for cq = 0.1, 0.3, 0.6, 1, 2, and 10 in ascending order from 
the green cq = line (green). 

x-bending solutions 

A similar analysis can be done for the ^/-bending solutions. The behavior of these solutions 
in the global Kruskal coordinates is easy to deduce from what we already know. The C2 = 
solution is the straight string which extends into the second asymptotic region. Likewise, 
all the < [c2 1 < 1 solutions extend to infinity in the second asymptotic region. The 
|ca| = 1 solutions asymptote to the horizon and extend to \y\ = oo. Finally, the \c%\ > 1 
solutions turn at a finite altitude above the horizon and extend back to infinity in the first 
asymptotic region. Since none of these solutions bend in the x°-direction, none of them hit 
the singularity. 

5.3 Euclidean signature 

We now do the same analysis in Euclidean signature. Since the OWL is at a fixed time, the 
Euclidean continuation of its one-point function is an insertion at fixed Euclidean time t E 
and should be interpreted as a finite temperature expectation value of the OWL operator. 

The Euclideanized metric is just (58) with dt 2 — > —dt 2 E 

d Jf = ^ldt 2 E + v?' 2 dx 2 + ^-du 2 + u^dQl (104) 

Near u — 1, this looks approximately like 
ds 2 1 

— w 3(u - l)dt 2 E + dx 2 + — -du 2 + dttl = dp 2 + p 2 d6 2 + dx 2 + dn 2 4 , (105) 

7^ " 3{u — 1) 

where in the second equality we changed coordinates u —> p := 2y/u~^l/V3 and t E -»• 9 ■ = 
3t E /2, since t E has period 47r/3 for the geometry to be smooth at u — 1. 

x°-bending solutions 

With the gauge and embedding 

x 4 = a\ u = a 2 , x° = t E {u), (106) 
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one finds that the action is given by 



From this, we find 
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if 



For c 7^ 0, solutions all reach the turning point at u 3 = 1 + Cq > 1 at finite (angle) t E , 
at which point the solution then goes back to large u at some new asymptotic value of t E - 
There is no value of Co for which the behavior changes qualitatively. At Co = 0, t E = constant 
down to u = 1. Near u = 1, this simply means that it is a straight line in the {p, 9} polar 
coordinates approaching the origin at a constant value of 9, say 9 = 0. Thus, it continues 
through the origin coming out at 9 = 7r. Figure 18 shows the value of A9, the difference 
in the asymptotic values of 9 = 3t E /2 for large u, as a function of c . Since all of these 
solutions turn and head back out to large u where there is no OWL operator insertion for 
them to end on, they should be discarded as solutions (or considered to have infinite area). 




Figure 18: AO as a function of c . The dashed line is (6v / 7rr[|])/(c ( 1 ) /3 r[i]), which is the asymptotic form 
of A6(c ). 



x-bending solutions 

As before, we choose the gauge and embedding x A = a 1 , u = a 2 and x 2 = y(u). Since the 
time coordinate does not come into play anywhere, the result is identical to the previously- 
discussed Lorentzian case, with the same qualitative behavior of the string worldsheets. 
Namely, for |c| < 1 the string goes through the u = 1 origin and comes out at A9 = it (and 
a different value of x as well). For |c| = 1 the string never makes it to u = 1 but instead 
asymptotes to u = 1 while extending infinitely along a x oc c direction. For |c| > 1 the string 
stays at 9 = and turns back up in u before reaching u = 1. 
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5.4 Interpretation of the curved worldsheet solutions 



The integration constants c M that appeared above are adjustable properties of the solutions, 
and should therefore reflect a change in the operator whose one-point function is being 
evaluated. Thus, instead of a unique operator OWL*(a; M ), as defined by Aharony and Ku- 
tasov, there must be an entire family of similar operators, OWL* c M ). What are these 
generalized OWL operators in holographic NJL? 

From the gravity side, it is not too hard to see what they are. Non-zero c M imply that the 
string worldsheet no longer ends perpendicularly on the Wilson line at u — u\. Equivalently, 
the Cfj, are proportional to non-zero fluxes of spacetime energy- momentum flowing into the 
string. This means that there are forces on the ends of the string. 

Since strings end on D-branes, in order for there to be a precise string dual of the 
Aharony-Kutasov OWL operator, it should be thought of as a contour on an appropriate 
probe D4- or D6-brane (at u = u\ and parallel to the N c color D4-branes) on which the string 
ends. In this case, the source of the forces on the string endpoints when c M 7^ is apparent; 
they are constant electromagnetic fields oc S^c^ on the probe D4- or D6-brane 3 . 

Some interesting things emerge from this picture. In the Lorentzian case, turning on c 
corresponds to an electric field on the probe D6-brane. |co| = 1 corresponds to the critical 
electric field beyond which the D6-brane no longer exists, and such an operator cannot be 
defined. We have seen that the corresponding |c | > 1 string worldsheets are those that hit 
the black hole curvature singularity. Likewise, turning on c corresponds to a magnetic field 
on the D6-brane. Here, there is no critical value but one expects that, for |c| > 1, there is 
a qualitative change in behavior as the strong magnetic field acts to "confine" the D6-brane 
to a smeared D4-brane; indeed, a qualitative change happens here — the strings go through 
to the second asymptotic region for \c\ < 1 but turn and go back out the first asymptotic 
region when \c\ > 1. 

However, this pleasing picture on the string side does not shed much light on the de- 
scription of the 7^ OWL operators in the effective five-dimensional Yang-Mills plus four- 
dimensional fermion gauge theory. Since the D6-brane electromagnetic fields are constant, 
the modification to the Aharony-Kutasov OWL operator should be distributed uniformly 
along the length of the five- dimensional Wilson line. It presumably corresponds to "dress- 
ing" this Wilson line with a density of gauge covariant insertions. One possibility, which 
breaks the four-dimensional Lorentz invariance in the correct way is 



0WL](^,O = i$(af,-e /2) V exp 



V2 



(iA 4 + c v F u4 + §)dx A 



-f. /2 



(109) 



In any case, whatever their precise form, their expectation value is a xSB order param- 
eter. We have found that, in the parallel D8- and D8-brane geometry, all the corresponding 



3 This situation is similar to that of time-like strings in the D3-D7 system, for which the Wilson line is 
the path of the string endpoint on the D7 brane, while the integration constants are interpreted as forces 
applied to the ends of the string by constant electromagnetic fields on the D7 brane. In the present case, we 
are essentially describing the same thing except with spacelike, instead of timclikc, Wilson lines. 
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string worldsheets have infinite area and thus lead to a vanishing order parameter, consistent 
with being in the xSR phase. 



6 More on generalized OWLs and curved worldsheets 

In this section, we consider more general worldsheets which are dual to the OWL operators 
in which the fermions are located at two different points in the ^-directions. We first 
consider the case in which the left and the right-handed fermions are separated only in the 
x°-direction, and then analyze the case in which they are separated in the x-directions. 



x°-bending solutions 

Consider the following gauge and embedding 

x 4 = a\ u = a 2 , x° = t(x 4 , u). (110) 

Following [1], this embedding should describe a worldsheet whose holographic dual is an 
OWL operator with chiral fermions separated in the x°-direction. The worldsheet has Neu- 
mann boundary conditions along the directions of D8 and D8-branes and Dirichlet boundary 
conditions for directions normal to them. With the above choice of gauge and embedding, 
the worldsheet action reads 

Sf= 2^J dudx 4 ^G uu G 44 + GuGuit')* + G tt G uu (iy, (111) 

where in this section we use ■ and ' to denote derivatives with respect to x 4 and u, respectively. 
For the worldsheet to be well-defined everywhere, one has to impose the condition 



1 

7 



2 det(M = G UU G U + G tt Gu(t') 2 + G tt G uu {if > 0. (112) 



The Lagrangian is invariant under a shift in x°, giving rise to the equation of motion 



GuG uu t \ I GuGat' 



i 



\J G uu Gu + GuGuit') 2 + G tt G uu {i) 2 ) \ a/ G uu Gaa + GttG^{t') 2 + G t tG uu {i) 2 



0. 



(113) 

Although we do not know a closed form solution to (113), its solution can be obtained 
numerically. 

However, instead of following this route, we will restrict our ansatz further in the hope 
of finding an analytic solution to the equation of motion. Choosing 

x° = kx 4 + t(u), (114) 
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with k being a constant, drastically simplifies the equation of motion. Indeed, substituting 
(114) into (113) results in an the following first integral of the equation of motion 

- GttGut ' = c, (115) 

V GuuGu + G tt G uu k 2 + G tt G 44(f) 2 

where c is a constant of integration (not to be confused with c defined in (18)). 

Note, however, that the simple ansatz of (114) does not solve the Neumann boundary 
conditions for the worldsheet. One can consider the contour of (114) as a good approximation 
to the contours where the complications due to the edges can be ignored. In that sense, the 
ansatz (114) merits further study. Solving (115) for t' results in 

t' 2 = c 2 U * - fcV ~ 1} (116) 

( M 3_1)2( C 2 + M 3_ 1 )' ^ ii0 J 

and substituting (116) into (112) yields 

1 , , u 3 -k 2 (u 3 -l) 

—deth a p = — — 

7 Z <r + u 



2 detv= o^: 3 _ 1 >°- ( 117 ) 



The positivity of the determinant of the induced metric (117), together with the reality 
of the solution (116), implies that c 2 > 0. Now, c = does not represent a physical solution 
for the following reason. At = k 2 /(k 2 — 1) the determinant changes sign, implying that 
the allowed range of u is [%,«»). Therefore, for c = to be a solution, the worldsheet 
must either have a turning point at some u t > so that it never reaches u* or, if it does 
not have turning points, then it must end at u = w*. Equation (116) implies that the only 
turning point is at the horizon u — 1, which is actually below u*. Thus, the solution must 
end at u = before it hits the would-be turning point. However, there is nothing for the 
worldsheet to end on as there are no D-branes there. Therefore, we discard the c = branch 
as a physical solution. 

For c 2 > 0, the story is the same as the c = case. In order to obtain a valid worldsheet, 
(117) indicates that again the worldsheet must end at u = [k 2 /{k 2 — l)) 1 ^ 3 . Since there is 
no D-brane there for the worldsheet to end on, we discard these solutions as well. So we 
conclude that, modulo a subtlety having to do with the boundary conditions, since there is 
no classical worldsheet of the type (114) the vev of the dual holographic operator must vanish 
in the xSR phase of the model. It is reassuring that the vev of this operator vanishes in 
the chirally-symmetric phase because we argued above that the dual operator is a chirally- 
charged OWL where the left-handed and the right-handed fermions are separated in the 
x°-direction. 



x-bending solutions 

For an OWL with the left and the right-handed fermions separated in, for example, the 
y-direction, we choose the following gauge and embedding for the holographic worldsheet 

x A = a 1 , u = a 2 , y = y(x A ,u). (118) 
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The worldsheet action takes the form 

,2 



Sv = 



7 



2ixct 



dudx \ / G 22 



G uu (l + {yy)+G 22 {y>y 



The Euclidean worldsheet must satisfy the condition 



G 



22 



> 0. 



G uu (l + (y) 2 )+G 22 (y') 2 
Since the integrand in (119) is independent of y, the equation of motion is 

/ V / 



(119) 



(120) 



G 22 G uu y 



\ 



G 22 G uu ( 1 + (y) 2 ) + G 22 (y') 



V 



(G 22 ) 2 y' 



V 



G 22 G uu (l + {y)A+G 22 {y'f 



0.(121) 



In order to find a closed-form solution, we relax (118) and choose a simpler ansatz 

u = a 2 , y = kx 4 + y(u). (122) 



x 4 = a\ 



Again, this simple ansatz does not satisfy the Neumann boundary conditions but is still 
worth studying because many contours can be approximated by such an ansatz as long as 
the edges of the worldsheet do not drastically change the shape of the worldsheet. 



Substituting (122) into (121) yields the first integral of motion 

(G 22 )V 



= c, 



G 



22 



G uu (l + k*) + G 22 (y>y 



(123) 



where c is a constant of integration. Solving (123) for y' and substituting the values of G 22 
and G uu gives 



y' 2 = (1 + k 2 )c 2 



{u 3 - l)(v?-c 2 )' 



(124) 



Except for a factor of 1 + k 2 , (124) is identical to (the square of) (99). Thus, our analysis in 
the previous section goes through here as well. In particular, when c = 0, we have a straight 
worldsheet which passes through the horizon and hits the second asymptotic boundary. Since 
there is no OWL insertion at the second asymptotic boundary, this worldsheet solution has 
nowhere to end and so one should discard this solution. When c 2 = 1, the worldsheet 
asymptotes to the horizon. The area of the worldsheet is 



27TOi j 1 



UA 



du 



u 



vfi — 1 ' 



(125) 



which is infinite (besides the usual UV-divergent piece) giving zero contribution to the one- 
point function of the dual OWL operator. When < c 2 < 1, the worldsheet dips into 
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the horizon and has a turning point at some radial position below the horizon from which 
it turns back up hitting the second asymptotic boundary of the spacetime. Since there is 
no OWL insertion on the second asymptotic boundary for the worldsheet to end one, one 
should discard this solution as well. Finally, when c 2 > 1, the worldsheet has a turning point 
at u t = | C2 1 2//3 (above the horizon) and turns back up to the boundary of the spacetime for 
which there is no OWL insertion to end on. This solution is non-physical too. Therefore, one 
reaches the conclusion that, in the xSR phase of holographic NJL, the one-point functions 
of the OWLs described above (those in which the fermions are separated in the ^-direction) 
all correctly vanish. 
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